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One component of a massive black hole binary (MBHB) might capture a small third body, and 
then a hierarchical, inclined triple system would be formed. With the post-Newtonian approximation 
including radiation reaction, we analyzed the evolution of the triple initially with small eccentricities. 
We found that an essentially new resonant relation could arise in the triple system. Here relativistic 
effects are crucial. Relativistic resonances, including the new one, stably work even for an outer 
MBHB of comparable masses, and significantly change the orbit of the inner small body. 
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Merger of two massive black holes (MBHs) is one of the most violent events in the universe. A huge amount of energy 
' is released in the form of gravitational waves (GWs). For a MBH binary (MBHB) with two masses ~ lO^-lO^M©, 
the Laser Interferometer Space Antenna [l[ can easily detect the waves from virtually anywhere in the universe. 
Meanwhile, GWs from an extreme-mass-ratio inspiral [EMRI, more specifically, a compact object {e.g. neutron star, 
white dwarf) orbiting around a MBH] would enable us to closely examine gravitational theories using a map of the 

■ spacetime around the MBH, encoded in the waves 0,3 . 
In a MBHB, either of MBHs might trap a small third object and form an EMRI {e.g. Such a compound 

O . triple system would be intriguing for cosmology and astrophysics. For example, with a tidal disruption of the small 
body, the MBHB merger could have precursive electromagnetic-wave signals that might allow us to identify the host 
galaxy of the MBHB and its redshift. Then we could observationally constrain the dark energy, using the luminosity 
distance estimated from the measurement of strong GWs emitted by the MBHB . 
^ ] In the past ~ 200 years, orbital resonances have been ubiquitously identified among planetary or satellite systems 
obeying Newtonian dynamics j8|. For example, the mutual stability of Pluto and Neptune is maintained by their 3:2 
orbital periods around the Sun. Therefore, one might expect that orbital resonances can be an effectual mechanism 
to form a stable three-body system including two MBHs. Indeed, except for detailed points, it was found that the 

^ ^ resonant relations similar to those known in planetary/satellite systems could set up relativistic triple systems evolving 

\^ . by emitting gravitational radiation (see e.g. Q for first order resonances such as 2:1 and also [l^] for co-orbital ones). 

■ As in the case of Newtonian systems, the mass ratio of a MBHB involved in these familiar and strong resonances 
] should be much smaller than unity (typically ^ O(10~^)) [ll|. However, from the viewpoint of gravitational wave 

■ astronomy, it is preferable that a MBHB of comparable masses can resonantly trap a small third body. 

] In this paper, we report that an essentially new resonant relation could arise in a compound EMRI/MBHB system 

■ and it would significantly change the orbit of the inner EMRI. Here, relativistic effects, and inclined, hierarchical 

■ orbital configuration of the triple are crucial. These would be naturally realized for an EMRI/MBHB system, and 
^ , the new resonance is among the two strong relativistic ones that work even for outer MBHBs of comparable masses, 

• ■ unlike the observed planetary or satellite systems. 

rS [ This paper is organized as follows; in §11 we describe our numerical method and summarize basic notations. In 
^ ■ §111 we show a typical orbital evolution under the new resonant relation. §IV is devoted to studies on the resonance. 
■ " " ' We first examine analytically how and why the resonance appears in our simulations, and then provide results from 
systematic numerical analyses. We also discuss issues related to gravitational wave measurements. In §V, we roughly 
evaluate the expected capture rate of a small third body for a MBHB. §VI is a summary of this paper. 



II. NUMERICAL METHOD 



In this paper, we numerically study evolution of a triple system composed by a MBHB (masses mo = (1 — q)M 
and mi — qM) and a small inner particle m2(^ M) orbiting around toq (see Fig.l). Here M is the total mass of 
the MBHB and q is its mass ratio. We assume that niQ and mi are comparable. The two objects mo and m2 can be 
regarded as an EMRI progressively perturbed by the outer body mi whose distance to mo shrinks due to gravitational 
radiation reaction. Below, we adopt the geometrical unit with G — c — M — 1, and do not include effects of spins. 

For equations of motion of the triple, we use the ADM Hamiltonian at 2.5 post-Newtonian (PN) order (see e.g. 
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[T^ [isj ) formally written by 

H = Hn+Hi+H2+ H2.5- (1) 

Here iJAr, i/i and are the Newtonian, the IPN and the 2PN terms respectively. H2.5 is the first dissipative term 
induced by gravitational radiation. The Hamiltonian H is originally given for the position variables (i = 1,2,3: 
suffix for the particles) and their conjugate momenta Pj. But, instead of Pj, we introduce new variables Si = p^/rrii 
and take appropriate partial derivatives 



I dH . _ 1 dH 

m, ds,, ' * rrii dr,. 



(2) 



in order to improve accuracy at numerically integrating systems with large mass ratios (the dot ' representing a 
time derivative). This prescription enables us to safely analyze the system even in the limit TO2 ^ 0. But the three 
bodies toq, mi and TO2 are handled equivalently in the post-Newtonian framework without introducing approximations 
associated with 7712 <Si mQ,mi. We numerically integrate the equations of motion using a fifth order Rungc-Kutta 
method with an adaptive step size control (see @ for details of the numerical method). 

We set up the initial conditions (denoted with the suffix "s") of the compound EMRI/MBHB system, as follows. 
First, for the outer MBHB, we put its semimajor axis at ais > 300 with the circular orbital velocity including the 
IPN correction. We can realize a small initial eccentricity eis = 0(10"''). Next, for the inner EMRI, we inject the 
small particle m2 at the distance a2s(<C ais) from the moving body mo, and set its relative velocity at the Newtonian 
circular velocity, to generate a small eccentricity (typically e2s = 0(10"^)). For simplicity, we mainly set the initial 
eccentricities eis ^ 1 and e2s ^ 1- Actually, for a single MBH, a majority of EMRIs might originate from dissolutions 
of stellar mass binaries by the MBH [13], and later become 62 ^ 1 due to radiation reaction. In §V we revisit issues 
related to the orbital parameters of the inner EMRIs. Meanwhile, assuming independent evolutions of EMRI/MBHB 
in earlier epoch, we randomly put their mutual orbital phase and inclination. 

To monitor the orbital elements, we define the coordinate distances between the three particles as 

di = dio = |ri - ro|, ^2 = ^20 = |r-2 - ro|, di2 = |ri - r2|. (3) 

For i = 1,2, the semimajor axis and eccentricity (including the forced ones) are calculated from the distance 
di through its consecutive maximum (0^(1 + e^)) and minimum (0^(1 — e^) = Vpf. pericenter distance). The MBHB 
evolves predominantly by its gravitational radiation alone, and we have ai ~ di with ei ^ eis(ai/ais)^^/^^ <^ 1 Jl5J. 
Since we perturbatively include the relativistic effects, we terminate our numerical integration when the distance dij 
between any pair i-j becomes less than 10 times rrii+mj (either CI: di < lOM or C2: d2 < lOmo or C3: di2 < lOmi). 

To define the angular variables of the inner and outer orbits, we introduce the Cartesian frame XYZ around the 
central body mo with fixed spatial directions (see Fig.l). Here the XY-p\a.ne is identical to the initial orbital plane of 
the MBHB. Because of the initial inclination of the EMRI, the orbital plane of the MBHB slightly precesses. However, 
given mo, mi ^ m2, the outer body mi virtually stays on the XY plane, and we use Ai for its angular position. 
Following standard conventions we also define ^2 and ru2 for the ascending node and the pericenter of m2. The 
inclination I2 is the angle between the angular momentum of m2 and the Z-axis. 

To clarify the relative position between mi and m2, we also introduce a complimentary frame XnYnZn that is 
corotating with mi. The rotating Xjj-axis is oriented from mo toward mi, and the Z/j-axis coincides with the original 
Z-axis. The coordinate values in the two frames are related by = Rz{—Xi)X . Here Rz{—Xi) represents the 
3x3 rotation matrix around the Z-axis with the angle — Ai (and so on below). 

For reference, we provide useful expressions for a binary with masses m and m', semimajor axis a, and the eccen- 
tricity e. The Kepler angular frequency n, the relativistic apsidal precession rate by the IPN term [l6| are given 

by 

^^'^ 3n(m + m!) 

, = — — . (4 

a(l — e^) 

Meanwhile, the orbital parameters a and e decay due to gravitational radiation reaction as (l5| 

(5) 
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FIG. 1: Configuration of the compound EMRI/MBHB system. Tlie XY Z-bame is defined around the central body mo, and 
its orientation is spatially fixed. The massive outer body mi virtually remains on the XK-plane (its initial orbital plane and 
used as our reference plane) with its angular position Ai. The inclined inner particle m2 intersects with the XF-plane at the 
ascending node specified by ^2. Its pericenter is described by the angle 7172 = fl2 + 0^2 {i^2'- the angle between the ascending 
node and the pericenter). A2, fli and zui can be defined similarly, but less important in our study. 



III. ORBITAL EVOLUTION 



In Fig. 2, we show the orbital elements of the inner EMRI in one of our runs. The three masses are mo = 0.6, 
nil = 0.4 and m2 = 10~^. We set the initial orbital parameters, ais = 350, = 37mo, I2S = 74.3°, eis ^ lO""* 
and e2s 3 x 10^'^. In Fig. 2, we use the outer semimajor axis ai as an effective time variable moving leftward from 
fli = 350 down to ^ 160. 

At first, the EMRI evolves almost independently on the distant massive body mi, and its orbital decay rate 02 is 
close to the analytical prediction ([SJ. Then, at ai ~ 290, the eccentricity 62 and the inclination I2 start to increase 
rapidly. The orbit of the EMRI becomes retrograde at ai ^ 240. Our integration is ended at ai ^ 160 by the condition 
C2 (02(1 - 62) ~ 6), when the time before the merger of the MBHB is ~ 5.4(M/10^MQ)yr. 

In Fig. 2d, we plot the combination of the angular parameters (modulo 27r) 

(fi = 3Ai — n72 — 2^2 (7) 



sampled at intervals. This key variable initially shows no structured pattern, but becomes localized around ip ^ +0 
at 02 ^ 300 (satisfying \w2\ ^ 1^2 1)- It is clear that the EMRI is now resonantly trapped by the MBHB. To the best 
knowledge of the author, this is a new resonant state never discussed in the literature. As we see later, relativistic 
effects and the inclined, hierarchical orbital configuration are the crucial elements to raise the resonance. We also find 
that, for H = (without the PN terms), the eccentricity 62 of this highly inclined system is quickly increased by 
the Kozai process 17i] (easily destroyed by the relativistic apsidal motion J^) and the run is ended shortly by the 
condition C2. 

The combination ip indicates that the pericenter of the inner particle 1112 has a simple geometrical relation to the 
position of the massive outer body mi. In Fig. 3, we plot the snapshots of the pericenter and the position of TO2 
around the epoch ai ^ 265. We use the corotating frame Xj^YrZ/}, and the pericenter lies almost on a distinct 
one-dimensional structure. Indeed, the shape of the curve is roughly given by the simple analytical expression 

(cos[3w] \ 
sin[3ii] J (8) 

parameterized hy u ^ Xi — Q2 ^ {'^2 — ^^2)/3. Here the rotation Rx{l2) represents the tilt due to the inclination and 
Rz{~u) is for the conversion to the corotating frame. The rapid evolution of 62 and I2 in Fig. 2 can be understood 
as accumulation of coherent interactions between EMRI-MBHB though the established resonance. 
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FIG. 2: Evolution of the inner EMRI, as a function of the decaying semimajor axis ai of the outer MBHB (leftward from 
ai = 350 down to 160). The panel (a) is for the inner semimajor axis 02. The inner eccentricity 62 (panel b) and the inclination 
I2 (panel c) show sudden changes around a2 ~ 290 where the inner EMRI is resonantly captured by the outer MBHB, as shown 
in panel (d). 
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FIG. 3: The positions (crosses) and pericenters (filled squares) of the inner particle m2 seen in the corotating XRYRZR-frame 
normalized by the outer distance ai. In this frame, the central body mo is at (0,0,0) and the outer body mi is at (1,0,0). 
The left panel is the projections of the points onto the XflVR-plane, and the right one onto the ^flYR-plane. These results are 
obtained from the same run as Fig. 2. The points are sampled 200 times between ai = 272 and 256, corresponding to ~3400 
orbital cycles of the outer MBHB. The pericenter of the EMRI is resonantly trapped by the MBHB, and stays nearly on the 
analytical curve ^ with = a2(l — 62). 



IV. RESONANT MODES 



A. Analytical Studies 

Now, with the aid of the disturbing function, we discuss why the resonant variable if appeared in our triple system. 
Roughly speaking, the disturbing function is a perturbative expansion of the effective gravitational potential between 
TOi and 7712 both orbiting around a central body mo, and given by a summation of terms proportional to 

cos jiAi + ki-cu^ + kQi^ (9) 

with integers ji, hi and U Q. The disturbing function is usually applied for systems with mg 3> mi, 7712 {e.g. planets 
around a star), but would prove to be quite useful to interpret our numerical results with mg ^ mi S> ?n2- 

Since the outer MBHB continuously has ei ^ 1 and its orbital plane is almost unchanged, its pericenter and 
ascending node would not be important for the present resonant capture. Therefore, in the disturbing function, we 
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analyze the terms Aj-^j^ki cos (jjj-^j^^ki with the phases 

<^Jli2fci = Ji^i + i2A2 + fcro2 + lft2- (10) 

Below, we assume ji > 0, as we are interested in explicit resonant couplings between EMRI-MBHB. 

We rely on the three basic properties generally valid for the individual terms of the disturbing function 0, [l^ ; 

(i) the sum rule; ji+j2 + k + l = from the rotational symmetry around the Z-axis, 

(ii) the scaling relation of the amplitude; Aj-^j^ki — O Cj'^'/j' (02/01)^ for 62, /2, (02/^1) ^ 1 with p > 2, 

(iii) the restriction of I to even numbers: I — 1v (v: integer) from the symmetry with respect to the XK-plane. 

A resonant state is identified by the condition (t>j-^j2ki — const. Here we analyze a more tractable form 4'j-^j2ki ^ 0, 
and evaluate the magnitudes (Ai, A2, n72, ^^2) and the adequate integers {ji, i2,k,l) for the triple system with the 
hierarchy 02 ^ ai and a small initial value e2s ^ 1- 

Because of the correspondence A^ ~ cx ^^"^ (see eq.Q) and the general relation A^ 3> Wi, (li valid also in weak 
field regime, our system satisfies A2 3> Ai, n72, ^2- Then, for j2 7^ 0, the condition (j)j-^j2ki ^ implies a very high-order 
resonance with \ji + ^2! ^ \ji\ ^ |j2|(ai/a2)'^^^ ^ 1- Instead, we limit our analysis only for the simpler cases with 
J2 = 0. 

Next we compare 1372 and f22, by separately evaluating the secular Newtonian effects from the distant body mi and 
the relativistic effects around the nearby one mg. The Newtonian contributions '!i7N2 and i}N2 become 0(712(02/01)'^) 
Q, while the IPN ones are 117112 ^ 3ri2(l — q)/a,2 (see eq.(|3])) and f2fl2 — [s^]- The total relation becomes 
W2 — wr2 ^ 1^22! ^ |^^7V2|, consistent with our numerical results. Thus, for our weak field system, the resonance 
fjij2ki — const [ji > 0) should be realized with the relation jiAi + fcri72 ^ (A; < 0). Here, a smaller \k\ is preferred 
from the scaling relation (ii), and we put k = —1. In contrast, a small inclination angle I/2I ^ 1 is not assumed, 
and we do not need to impose a strong requirement on I = 2^ at present. Then, the resonance variable valid for our 
system is written as 

e.,j, = {2u + l)Xi ~ m2 ~ 2^2. (11) 

Below, we attach ";" before the mode-number suffixes 1/ to distinguish them from the labels i for the particles. 

As seen so far, the relativistic effects are crucial to make 9-i, ~ by increasing the apsidal precession rate W2. In 
addition, the overall dissipative evolution is due to the 2.5PN radiation reaction force. At the resonance, the PN order 
parameter (1 — q)/a2 = 0(^72/77.2) of the EMRI is comparable to the orbital hierarchy (02/01)'^/^ ^ 711/712. Thus, our 
sequence O-i, contains only the outer position Ai without the inner one A2, remarkably different from the standard 
mean motion resonances with jij2 7^ 0. Here, the unusually large outer mass 7711 (^ mo) would enhance the resonant 
coupling. 

The observed combination ip in eq.(I7|) is properly reproduced as (p = 9-i. Interestingly, for v — 0, the variable 9-0 
coincides with the j = 1 inner Lindblad resonance. The Lindblad resonances [here parameterized by j'Ai — (j — 1)A2— 'Ci72] 
play fundamental roles in the dynamics of coplanar disks [l^l whose relativistic effects are recently discussed in [21j] . 

From a basic relation of the disturbing function {e.g. eq.(12) in [lOj), we get p — max{3, 2i/ + 1} for the relation 
(ii) with the variables 9-^. This shows weaker coupling for v > 2. Actually, even for \k\ — 2, we have p < 3 only for 
(j) = 2Ai — 2tx72 that completely degenerates with 9-o = Ai — Tx72. 

In the C ring of Saturn, there is a ringlet structure whose longitude of the pericenter tu is in a resonant relation 
with the angular position At of Titan, the largest satellite of Saturn. The resonant variable is given by At — and 
the apsidal precession w of the ringlet is mainly caused by the large multipole moments of Saturn {e.g. its quadrupole 
moment J20 = 0.016) [l^l, instead of relativistic corrections. 

Now we can predict when the inner EMRI with 62 ~ is resonantly captured by the outer MBHB. From the relation 
(27/ -|- l)Ai — -072 — (or equivalently 27/ + 1 ~ Ti72R/'ni), the critical inner semimajor axis 7;j/ is given as a function of 
ai by 



32/5(27/-Hl)-'/'(l-<7)3/5af^ (12) 



For u — Q and the test particle limit q = nii — > 0, this expression coincides with eq.(130) in [21] for the j = I inner 

Lindblad resonance. Using the results in Fig. 2, we also confirmed that the ratio -072/711 cx aj^^a2 '^^^/(l — 62) is nearly 
constant during the trapping. 

B. Numerical Studies 

Next we systematically analyze a series of numerical simulations. We take the mass parameters; q = 0.4, ?Ti2 = 10~^ 
and the initial conditions; eis < 10~*, e2s < 10^^, ais = 350. We prepared totally 11 runs from various initial inner 
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FIG. 4: Evolution of the triple systems in the aia2-plane. The open squares are the initial conditions with ais — 350 and 
a2s = 20-30. The circles represent the points when 62 exceeded 0.1 (indicating resonant capture), and the triangles are the 
termination points of the runs (all by the condition C2) . The two solid curves are the analytical predictions (|12|l for the onsets 
of the resonant capture for 1/ = 1 and 0. 

separations a2s G [33.3mo, 50rno] with random initial inclinations cos /2s G [0,1]- For these sets (ais,a2s), the infall 
rate is larger for the MBHB (namely d{a2/ai)dt > 0). From eq.([5|), this catch- up condition is written as 02 > f^cOi 
with rjc = {mom2/miAiy/'^ . When we reverse the time, each non-resonant triple moves on the curve 



and asymptotically approaches to the line 02 = ?7cii- 

In Fig. 4, we provide the time evolutions of the runs. We obtained similar results without the 2PN term. The run 
from a2s — 22 is what was already shown in Fig. 2. In Fig. 4 we added the analytical predictions 02 = "f-u (y — 0, 1) for 
the onset of the resonant capture. They reasonably agree with the numerical results. Since the initial outer distance 
a\s = 350 is not sufficiently large for a2s ^ 25, the corresponding EMRIs are captured by the i/ = mode. We also 
examined dependence of the resonant captures on the mutual inclination of the two orbits. It was found that slightly 
inclined EMRIs (e.g. /2s = 0.14 and 0) could pass through the v = \ mode and first reacted to the v — ^ mode, in 
accord with the scaling relation Aj^j^ki oc {hY" ■ 

We briefly describe other interesting results. First, to realize a capture, the resonant curve 02 — 7;i/ {v — 1,0) in 
Fig. 4 should be encountered in the direction v < Q. For > 0, both 62 and I2 show gaps at the resonant crossings, 
but the capture was unsuccessful. Using eg. (1121) in the form {2v + 1) cx a^''^a2 we can assign contour levels v on 
the (01,02) plane. Then, from eq.dS]), we obtain 



with K = (5/3)1/''. Therefore, we have v <0 (required for captures) at 02 > nijcai. In Fig. 5, we provide a schematic 
illustration for the resonant capture. We denote the intersection of the line 02 = nrjcai with the resonant curve 
02 = 7;j/ by (a;i;i/, a;2;j/). For a capture by the v mode, a triple should cross its resonant curve at Oi < Xi-^ (the solid 
part in Fig. 5). Furthermore, as understood from the flows of triples in the (oi, 02)-plane, an EMRI capturable either 
by = or 1 must exist in k~^X2;o < 02 < X2-o (shown by the double line in Fig. 5) at the critical epoch oi = xi-q. 

Secondly, for the parameters shown in Fig. 2, we examined systems with larger 625. For 62 0.11, the system 
transversed the v = 1 curve without a capture but again showed gaps of 62 and J2. In the same manner, we analyzed 
the crossing of the v = mode around oi ~ 200 and 02 ~ 29 with various eccentricities 62, and found an upper limit 
62 ~ 0.35 for yielding captures. Similar behaviours (e2- or sign(t')-dependence of the captures, observed gaps of 62 at 
the resonant crossings without captures) are found for standard mean motion resonances and well explained by the 
separatrix structure in an effective phase space (see Figs. 3-6 in ^9|] and also j2^). 

By dropping the time consuming 2PN term, we also performed runs from ois ^ 600 and 02^ ^ 30 (corresponding to 
I' > 3 in Fig. 4). But the captures by the v > 1 modes did not occur, as anticipated from the arguments on the power 
p in the relation (ii). Furthermore, by artificially multiply a large factor to the dissipative term H2.5, we evolved 
widely separated systems (oi ^ 2000 and 02 ^ 7;o ^ 100 for = 0) at accelerated rates, and confirmed captures by 
the V = mode. 




iVc'^i)^ — const 



(13) 




(14) 
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FIG. 5: Resonant capture of the inner EMRI by the outer MBHB. We show a schematic iUustration for the v — mode with 
a smaU eccentricity 62 — 0. The thick (sohd and dashed) curve is the critical semimajor axis 02 — j-o of the EMRI for the 
resonant condition (|12|l . The solid line 02 — ri^ai is the catch-up line where we have a^^dai/dt = a2^da2/dt. The dotted line 
a2 = Kr^cCLi with k = (5/3)^'''' divides the signature of the resonance approach sign(!>), as shown by the three arrows. We put 
the intersection of the two curves by {xi-o,X2-o)- To realize a capture, the resonant curve 02 ~ 7;0 should be crossed in the 
direction i> < 0, corresponding to the upper region of the dotted line (the solid part of the curve 02 = 7;o). Therefore, at the 
critical outer distance ai = a;i;o, the capturable inner EMRI must has a semimajor axis k~'^X2;0 < 02 < X2;0 shown by the 
vertical double line. 



C. Gravitational Wave Measurements 

Here, we discuss GW detection for a trapped EMRI. In Fig. 2, the final infall rate fp2 = d[a2{l — e2)]/dt is ^ 70 
times larger than that of a corresponding isolated EMRI (see eqs.(l5]) and Considering the effective time duration 
of the GW signals, the detectable distance for the trapped one becomes ^ I/VTO times smaller. This is partly due to 
the small inner mass m2 = 10^^, and, indeed, the results in Fig. 2 are similar to those for a test particle with m2 = 0. 
We thus analyzed systems with larger m2 from (ai^, a2s) = (305, 21.5) now down to — 6mo. The resonant trapping 
was successful up to m2 ~ 10~^ for which the trapping ended at rp2 ~ 8mo with the final rate fp2 close to the isolated 
EMRI. But the corrections higher than 2.5PN could become important here. 

To detect a trapped EMRI, a simple data analysis worth trying is a search triggered by strong GW of a merging 
MBHB that provides e.g. mi, 7712, the sky location of the system and also a\ as a function of time. Here the primary 
orbital parameters of the EMRI would be its semimajor axis 02, eccentricity 62 and mutual inclination I2. These 
should be searched by using templates. The trapping condition a-^^ oc a^^'^iX — e|) might become useful to narrow 
the parameter space of 02 and 62 to be surveyed. 

Now we go back to the specific system shown in Fig. 2, as an example, and discuss its final phase around oi ~ 160 
and 02 ^ 19. The orbital periods of the EMRI and the MBHB are ~ 700 and ~ 13000 respectively. Since the 
EMRI has a large eccentricity, it emits relatively strong pulse-like GWs around its pericenter. These GWs would be 
an interesting observational target. The characteristic duration of each pulse is (?'p2/™o)^^^ ~ 20 and the interval 
between the adjacent pulses is approximately the orbital period of the EMRI. Therefore, the semimajor axis 02 is the 
critical parameter for matching the GW signals. 

For a standard mean motion resonance, the inner and outer semimajor axes are related through a simple linear 
relation between the two orbital periods. For example, we have 02/01 ~ + 1)]^^'^ for the first order resonance 
with the librating variable JA2 — (.7 + l)Ai ^ const (j > 1: an integer). This kind of relation would be quite 
helpful to estimate evolution of the inner axis 02 from the observed GW signals of the merged MBHB. However, 
our hierarchical resonances 6-1, do not explicitly depend on the angular position A2 of the trapped EMRI, and we 
cannot directly deduce the key parameter 02 for the EMRI, as a function of ai. Therefore, for an EMRI trapped in 
our resonances, the identification of its GW signature could become more demanding, compared with standard mean 
motion resonances. 
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V. CAPTURE RATE 

In this section, we roughly discuss probabihty of a MBHB (2 x 10^ Af© + 10*^ Af©) trapping a lOM© BH at the ly ^ 
and 1 modes. With respect to these mass parameters, the critical point in Fig. 5 is given by 

(a;i;o,a;2;o) = (2040,118). (15) 

The formation scenario of an EMRI around a single MBH (not around a MBHB) often studied in the literature 
is the capture of a compact object around its close approach to the MBH by emission of gravitational radiation 
(hereafter two-body capture process) [23 - l26j . However, the formed EMRI would be highly eccentric, and the residual 
eccentricity has a typical value 62 ~ 0.5 even when entering into the LISA band, corresponding to 02 = 0(10) 25]. 
Therefore, considering the thresholds for the inner eccentricity 62 discussed in §IV.B, the resonant capture by an outer 
MBHB would not be promising for the EMRIs formed by the two-body capture process. 

Miller et al. p3l discussed formation of an EMRI through a tidal dissolution of a compact stellar binary by an MBH 
(hereafter binary dissolution process). One component of the binary is expelled from the system as a high- velocity 
star [i^] and another one becomes bound to the MBH. The resultant EMRI could have a larger initial pericenter 
distance (thus a larger cross section for its formation), compared with the two-body capture process, and could has 
comparable or larger population. In addition, for the binary dissolution process, the residual eccentricity 62 becomes 
small 62 ^ in the LISA band, due to a long-term dissipative evolution of the orbit with emitting gravitational 
radiation. 

Here, following [l3], we evaluate the orbital parameters of the inner EMRI formed by the binary dissolution process, 
simply assuming that the EMRI is affected by the secondary MBH solely through the v = and 1 resonances. 

We first consider only the evolution of the EMRI without taking into account the secondary MBH mi. A compact 
stellar binary with an orbital separation abin and mass ~ 7712 is dissolved by the primary MBH mo at the distance 

/o \ 1/3 / \ 1/3 / \ -1/3 

d,,,,^^ ^-!!^) a,,„ ^ 8.8 f AU, (16) 

where we put the fiducial parameter abm=0.1AU given in |14|. We regard d2,tide as the initial pericenter distance of 
the EMRI and also use the typical value 62 ^ 0.98 quoted in [lj| for the initial eccentricity. 

Then, with eqs.© and ([5]), the merger time of the EMRI is given by ~ 5 x lO^yr (less than the age of the universe). 
At a2 = X2-Q = 118, we have a significantly reduced eccentricity 62 ^ 0.17 and the remaining time 1.1 x lO^yr (oc a|). 
When the outer MBH inspirals down to the critical separation oi = xi-o = 2040, the distribution of the inner EMRIs 
at 02 < X2;o would be in a steady state, from our assumptions. We can characterize the distribution of the EMRIs 
by their infall rate at i?i„ = 5 x 10~^yr~-'^ [24|. This rate i?m was originally given for the two-body capture process, 
but we use it for the binary dissolution process, following the arguments in 14]. Finally, the probability of a MBHB 
merger with a trapped lOM© BH can be evaluated as 

P = (5 X 10"V"^)(1-1 X 10V)(1 - 3/5) = 0.022, (17) 

corresponding to the expected number of EMRIs on the double line in Fig. 5. 

Actually, the above probability P contains the contribution of the EMRIs that would be captured by the 1/ = 1 mode. 
These EMRIs pass close to the lower bound 02 — K~^ai-o — 103.9 on the double line. Since we have X2-m oc l/(2z^-|- 1), 
the branching ratio of these modes is given by 

1 - 3""^ ; S"'' 80 : 1, (18) 

and dominated by the v = Q mode. For inner EMRIs with white dwarfs or neutron stars of m2 ^ 0{1Mq), the merger 
times become larger than the age of the universe for the previous input parameters ahm — O.lAU and 62 ~ 0.98. But, 
if we can assume small evolved eccentricities 62 ^ 0.2 at 02 = X2;o and the steady-state distributions normalized by 
5 X 10^* yr^^, the probabilities become P = 0(1) for these EMRIs. 
Thus far, we have considered only an isolated three-body system in a simplified manner. In reality, a pre-existing 
EMRI might be destroyed by other stars {e.g. scattered from the secondary MBH). But, at the same time, the 
secondary MBH could highly enhance the capture rate of EMRIs Rim at least, for those formed by the two-body 
capture process [4]. This is an interesting possibility and would be worth examined in the context of the binary 
dissolution process. 

In ongoing or planned wide-filed surveys for transient electromagnetic waves (see [28j for recent results), we might 
detect a tidal disruption event resonantly driven by a merging MBHB. For such a event, the signature of the orbital 
period of the MBHB might be found in the temporal structure of the emitted electromagnetic waves (see [2^ for 
related discussions). Here, an MBH more massive than ^ IQ^Mq can directly swallow a white dwarf without a tidal 
disruption, and we need to consider a main sequence star in this mass regime. 
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VI. SUMMARY 

Orbital resonances might be a potential mechanism to append a small body to an inspiraling MBHB, and to main- 
tain the compound EMRI/MBHB system. Based on the post-Newtonian approximation, we numerically examined 
such hierarchical three body systems evolved by emitting gravitational radiation. For mutually inclined orbital con- 
figurations, we found a new resonant state with the librating variable 6-i = 3Ai — W2 — 202- This resonant state, 
together with another state 9.q = Ai — n72, works relatively strongly for the triple systems. Here the relativistic 
apsidal precession is essentially important, and the post- Newtonian parameter at the capture becomes comparable 
to the hierarchy of the two orbits (a2/ai)^/^. In contrast to standard mean motion resonances known among plane- 
tary/satellite systems, these relativistic resonances can prevail even for outer MBHBs of comparable masses. During 
these resonances, the eccentricity and inclination of the inner EMRI increase almost monotonically, and its pericenter 
distance coiild go down blow ^ 5 Schwarzschild radii. 

In order to realize a capture into the resonances, an inner EMRI cannot have a large eccentricity. Therefore, EMRIs 
formed through the familiar two-body capture process would be difficult to be involved in the resonances. In contrast, 
the binary dissolution process can produce mildly eccentric EMRIs at the spatial scale in interest, and these EMRIs 
might be resonantly trapped by inspiraling outer MBHBs. 

The author would like to thank Takahiro Tanaka, Xian Chen and anonymous referees for helpful comments. This 
work was supported by JSPS grant 2074015. 
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